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ABSTRACT 

A self-similar solution for time evolution of isothermal, self-gravitating 
viscous disks is found under the condition that o' = a{H/r) is constant in space 
(where a is the viscosity parameter and H/r is the ratio of a half-thickness to 
radius of the disk). This solution describes a homologous collapse of a disk via 
self-gravity and viscosity. The disk structure and evolution is distinct in the 
inner and outer parts. There is a constant mass inflow in the outer portions so 
that the disk has flat rotation velocity, constant accretion velocity, and surface 
density decreasing outward as S oc r~^. In the inner portions, in contrast, mass 
is accumulated near the center owing to the boundary condition of no radial 
velocity at the origin, thereby a strong central concentration being produced; 
surface density varies as S oc Moreover, the transition radius separating 

the inner and outer portions increases linearly with time. The consequence of 
such a high condensation is briefly discussed in the context of formation of a 
quasar black hole. 
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1. INTRODUCTION 


Quasars (QSOs) are the most powerful objects that have ever existed in the universe. 
The emergence of quasars at high-redshifts, ^ 5, is thus crucial when considering the 
formation of astrophysical objects, notably of galaxies. The view is widely accepted that 
QSO phenomena result from mass accretion onto supermassive black holes. However, the 
formation process of seed black holes at high redshifts is not well understood at the present. 
There are two distinct lines of thoughts concerning this issue. One is based on considering a 
formation of a proto-quasar supermassive black hole after the formation of a host galaxy as 
the consequence of stellar mass loss and star encounters at the nucleus of the galaxy (Rees 
1984). The other rather assumes a galaxy-independent population of massive black holes 
(Carr, Bond, & Arnet 1984; Loeb 1993; Fukugita & Turner 1996). Under the latter picture, 
a question is how quasar black holes formed at high redshifts, > 5 — 10. 

Suppose a high density fluctuation with a mass scale of ~ IO^Mq began to collapse at 
high redshifts oi z ^ 10. Such a cloud acquires angular momentum through tidal torque 
in the course of a gravitational collapse. Resultantly formed a rotationally supported, 
self-gravitating disk. For a typical spin parameter, the angular momentum barrier is 
by roughly seven orders of magnitude larger than the Schwarzschild radius (Loeb 1993; 
Eisenstein & Loeb 1995). The problem is then how to get rid of angular momentum from 
the cloud so as to form a black hole. Radiation drag via the cosmic background radiation 
seems to have been at work a.t z > 100, but is effective only when the cloud is optically 
thin (Loeb 1993; Umemura, Loeb, & Turner 1993; Tsuribe & Umemura 1996). Afterwards, 
angular momentum in the cloud could be redistributed via gravitational torque rising from 
nonaxisymmetric perturbations (Paczyhski 1978; Lin & Pringle 1987; Papaloizou & Lin 
1989) and/or turbulent shear viscosity which could be associated with magnetic helds 
(Shakura & Sunyaev 1973). It is thus worth investigating how a self-gravitating, viscous 
disk evolves in the context of black-hole formation. Furthermore, this kind of study is of 
great importance, of course, when one investigates physics of galaxy and star formation. 

The basic equilibrium structure of accretion disks are now well understood, as long 
as we believe the standard model based on the a-viscosity prescription (Shakura & 
Sunyaev 1973). Nevertheless, it is not easy to follow its dynamical evolution, mainly 
because the basic equations for the disks are highly nonlinear, especially when the disk 
is self-gravitating (e.g. Paczyhski 1978; Fukue & Sakamoto 1992). To follow nonlinear 
evolution of dynamically evolving systems, in general, the technique of self-similar analyses 
is sometimes useful. Several classes of self-similar disk solutions were known previously 
(Pringle 1974; Filipov 1984),but all of them considered a disk in a hxed, external potential. 
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We are now concerned with dynamical evolution of a self-gravitating disk in a 
time-evolving, self-consistent potential. As far as steady, nonviscous rotating disks are 
concerned, there are plenty of works so far done. Mestel (1963) was the first to find a simple 
disk solution, in which physical quantities are integrated vertically with respect to the disk 
equatorial plane. Hayashi, Narita, & Miyama (1982) found two-dimensional, isothermal 
disk solutions with hnite temperature (see also Toomre 1982 for stellar systems). Numerical 
steady solutions are calculated by several groups (Hachisu, Eriguchi, & Nomoto 1986; Bodo 
& Curir 1992; Hashimoto, Eriguchi, & Muller 1995). Recently, we have found a simple 
analytical solution for a steady, self-gravitating, isothermal disk (Mineshige & Umemura 
1996, hereafter Paper I) as an extension of Mestel (1963) disk. However, little study has 
been done concerning dynamical evolution of self-gravitating, viscous disks. 

We, in the present study, seek for a time-dependent, self-similar solution for a 
gravitational collapse of a rotation-supported, self-gravitating viscous disk. When a disk 
is sufficiently cool, gravitational instability will occur (Toomre 1964), providing a source 
of disk viscosity (Paczyhski 1978; Lin & Pringle 1987) or causing disk fragmentation 
(e.g., Bodenheimer, Tohline, & Black 1980). Several authors thus mainly discussed the 
consequence of gravitational instability in the context of fueling to active galactic nuclei 
(e.g. Shore & White 1982; Shlosman & Begelman 1987; Shlosman, Frank & Begelman 
1989), or (multiple) star formation (see Boss 1986; Myhill & Kaula 1992). We here adopt 
a rather distinct approach; we, in the present study, try to hud an analytical solution for 
a collapse of rotating, viscous disks, putting aside for the moment the stability argument. 
It might be noted in this context that Shu (1977) found the self-similar solution for a 
gravitational collapse of an isothermal sphere. Saigo & Hanawa (1996) discussed the effects 
of rotation. We extend these works so as to incorporate the effects of mass accretion via 
viscosity. We derive self-similar solutions in section 2, and then discuss the formation of a 
primordial quasar black hole in section 3. 


2. SELF SIMILAR, SELF-GRAVITATING DISK 
2.1. Basic Equations for Self-Similar Variables 


We start with the time-dependent version of the height-averaged equations for 
isothermal accretion disks (cf. Houma, Matsumoto & Kato 1991; Narayan & Yi 1994); 


(9S 

'm 


1 ^ 

r dr 


(rSW) = 0, 


( 1 ) 
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dt ^ dr 


_ GMr ^ 

p dr r 


^ 9{rV^) ^^_d_ . 3 ^ 

dt '' dr rT, (9r y dr J 

Here, S = 2pH is surface density, H is half-thickness of the disk, hi = V^/i?, Cg is sound 
velocity (which is constant by assumption), Mr is the mass of a disk within a radius r, and 
we approximated a potential to be ~ —GMr/r. This is a good approximation if S(r) profile 
is steeper than 1/r (see Appendix). We prescribe kinematic viscosity as 

V = aCsH = «(—)cgr, (4) 

r 

with a being viscosity parameter, because we will hnd later that self-similar solutions exist 
if a' = a{H/r) is constant in space. From now on, therefore, we assume a' (instead of a) to 
be constant. For vertically self-gravitating disks, H is determined as 


(47rGp)V2 27rG'S’ 2H Cg^ 

To proceed, it is convenient to rewrite mass conservation (1) using Mr{r,t)] 


dMr , dMr _ ^ 
dt ^ dr ~ ’ 


= 27rrS. 


Now, we introduce the following self-similar variables (Shu 1977); 

^ = Mr{r,t) = ^m{x), 

T/ / ^ V ^ 1 rV^{r,t) 

V^{r,t) = Csv{x), j{x) =xv = ——^ . 

Cg^ t 


Note that derivatives are transformed into 


d X d d d X d , , 

dt t dx dt' ’ dr r dx' 

for the transformation, (r, t) — {x,t' = t). Since all the time derivatives with respect to t' 
disappear if we use self-similar variables (Eq. 7), we hereafter write d/dx instead of d/dx. 

Equation (6) now becomes 

dm dm 

m + (u — x)—— = 0, and —— = xa, 
dx dx 


( 9 ) 
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yielding a simple relation between m, a and u] m = xa{x — u). With this being kept in 
mind, equations (1) - (3) can be modified as 


,1 da du u — X 
u- x)-— + — 

a dx dx X 


0 , 


( 10 ) 


2 da 
a dx 


+ {u 


x) 


du u — X 
dx X 
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0 , 


( 11 ) 


dj ,ld 
j + {u-x}— = a — — 
dx ax dx 


( 9 -^ 

ax —2j + X— 

V dx, 


( 12 ) 


2.2. Solution in a Slow Accretion Limit 


In the limit of slow accretion {v 3> l,a3> 1 ,|m| 1), equation (11) gives 


V = a^^^(x — j = a^^‘^x{x — 


(13) 


leading to 


dlnj 

din a; 


1 + 


1 1 


2x — u 


X — 


du '' 
dlnx, 


+ 


1 din a 


2 dlnx 


Note that from equation (10) we derive 


d In a 1 du 

dlnx x — udlnx 


(14) 


(15) 


from equation (12). Inserting equation (15) into equation (14), we have 


din j 
dlnx 


1 1 X 

2 2 X — M 


2x — u 
2(x — u) 


After some algebra, we obtain 


u , 1 d 

— = —a -— 

2x axj dx 


axj- 


2x — 3u 


X — u 


(16) 


( 17 ) 


With a help of the expressions for j (Eq. 13) and a (Eq. 15), we finally derive an ordinary 
differential equation for m(x): 


du 4x^ — 6ux + 3u^ 1 u{x — u^ 

dx 2(x — 3u)x a' (x — 3u)x 


( 18 ) 
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Equation (18) can easily be integrated numerically for an appropriate boundary condition; 
n = 0 at a; = 0 if we assume no central object (such as a black hole). Once u = u{x) is 
obtained, we can derive a = a(x) by integrating equation (15) for a given ag = cr(x = 1). 
The results of the integration are displayed in hgure 1 for different values of a' = 10“^, 10“^ 
and 10“^. The azimuthal velocity is derived from equation (13). 

Note that each physical quantity is a rather smooth function of x. We generally hnd 
du/dx < 0; that is, u{x) is a monotonically decreasing function of a; = r/cgf. Furthermore, 
physical quantities, such as u and a, are power-law functions of radius in the limits of 
x^ a' and x a'. 


In the limit of large a;(3> a'), mass accretion is induced by viscosity. Two terms on the 
right-hand side of equation (18) are balanced with each other (while du/dx = 0). We hnd 

u ~ —2a/ a ~ (Tox~^, v ~ m ~ 2a'ao, (19) 

where m(= —xau) corresponds to a mass-how rate. The radial dependences of physical 
quantities at large x are the same as those of the stationary, self-similar solution of a 
self-gravitating viscous disk (Paper I). However, we hnd W ~ —2aCs{H/r) in the current 
time-dependent solution, whereas T) = —OiCs{H/r) in the steady solution. This indicates 
that accretion velocity is doubled when we consider the ehects of continuously growing 
central mass (see discussion in Paper I). 


In the limit of small x ^ a' the hrst term dominates over the second on the right-hand 
side of equation (18), 


u —2x 1 — 


9 X 


11a' 

- 1/3 / 2 ; 


a 


(To f X 


- 5/3 


1 + 


8 X 


11a' 

1/3 




( 20 ) 


Note that u (and therefore W) is not proportional to a', indicating that mass-inhow is not 
controlled by viscosity, but is regulated by the inner boundary condition of W = 0 at r = 0. 
Mass is thus being accumulated continuously near the origin. 


To sum up, the disk structure and evolution is distinct in the inner and outer parts. 
The transition radius (rt^) separating these two parts increases linearly with time, because 
rtr ~ a'cgt oc t for a hxed a' (Eq. 7). We thus assume rtr = 0 initially; in other words, we 
consider the later evolution of the disk with S oc r~^ everywhere. (This is the situation 
postulated in Paper I.) As matter accretes towards the center, S prohle changes from inside. 

Now we recover physical variables from self-similar ones using equation (7): we obtain 

Vr ^ -2a'cs, S So ("—) , Vp ~ (27rGSoro)^/^, M ~ dWcsroSo, (21) 

\ro/ 
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at large r/t (3> a'cs), and 










14 ~ Cs 


/ r \ ^ 

\~o) \To) ’ 


M ^ dTrroSoCs 




( 22 ) 


at small r/t (<^ a'cs). Here, M = —27rrS14- is a dimensional mass-flow rate and we 
approximated M^. ^ /'’27rroSo(ir = 27rSorQ in equation (21). The units are 


ro = 1.0 Tpc pc, Cs ~ 10 ®'°T 4 '^^ cm s 



IqS.o 


pc 


ml/2 

-^4 


yr, 


Mo = dTra'csroSo ~ yr-\ (23) 

rpc 

for temperature of 10 "^T 4 K, mass of rsj 10^MqMq, respectively. The unit for S is chosen 
so as to give M = /27rE(r) rdr for the initial state, in which S = Egro/r; For such 
normalizations, a normalization constant of cT(a;) is 


En = 


M 

27rrn 


1 c iVfg 9 

g cm-2. 


r. 


pc 


_ 27rGfom 

(Jo =-z^o ~ 

c« 


10 


1.71 


Mo 

rpcT 4 ■ 


(24) 


Note that (Jq represents the ratio of disk radius to height at a; = 1 (see Eq. 7), or the initial 
ratio of gravitational energy to thermal energy of the disk, jc^ (Eq. 21). The model 
parameters of the self-similar solutions are o;', Cg (or temperature), and (Tq. 

Figure 2 plots the time evolution of a self-gravitating disk. Clearly, there are two 
regimes as mentioned previously (cf. Fig. 1). The radius separating the outer and inner 
parts is increasing linearly with time. If we follow a disk evolution at a fixed r, hence, we 
see that is initially constant and then decreases at t > r/a'cs- Accordingly, mass inflow 
rate also decreases with the time, causing a rapid growth of E and M^. Note that since 
H/r ~ {xa) ^ (Eq. 7), H/r is constant at large r/t, while it rapidly decreases inward; 
H/r oc (r/t)^'^. The thin disk and slow accretion approximations are even better in the 
inner portions at later times, although a may exceed unity at a; a'. This means, the 

present solution does not give a good representation of the disk structure at r/t q;'cs 
( discussed later). 


3. DISCUSSION 



3.1. Summary of the Self-Similar Solution 


We have derived a self-similar solution for time evolution of an isothermal, self- 
gravitating, viscous disk in the slow accretion limit. Disk structure changes from the inner 
to outer parts. For example, surface density is scaled as in the inner, while it is r~^ 
in the outer. This interface gradually moves outward in proportion to t. In this solution 
density increases monotonically with the time at the center. The mass prohle near the 
center is 

pr / T \ t T 

Mr{r) = J 27rS(r) rdr ~ 3 x 10® ^— j j —^Mq. (25) 

[Although this yields a diverging Mr, the increase of Mr should be terminated in a realistic 
situation, when the outer disk is depleted with gas.] 

As claimed hrst by Mestel (1963) and also by Paper I, the thin-disk approximation 
breaks down at radii comparable with the thickness. In fact, the present solution gives 
diverging D and a as a; approaching 0, which suggests that the solution does not represent 
physical situation at a; 1. Moreover, since we assume steady mass input towards the 
center, the central mass condensation increases at any time. Once a central object forms 
from a central mass condensation, gravity is dominated by this object at sufficiently small 
radii, where we may adopt a solution for a point-mass potential. 

Realistically, there may be two or three zones in a disk. Before forming an object, a 
self-gravitating disk has two zones (as mentioned in previous section). After the formation 
of a central object, in contrast there are three zones; the innermost region is dominated by 
a point-mass potential and the other two zones are dominated by self-gravity of the disk. 
Since M > 0, the mass of the central object is continuously increasing with time. The 
transition radius between the innermost to the inner region again increases linearly with 
time (Paper I). 

Self-similar solutions assume that heating and cooling rates have the same radial 
dependence (see Eq. 4 in Paper I). A flat temperature distribution is the result of this 
assumption. This is a reasonable approximation at least in the outer regions: when we 
balance viscous heating and radiative cooling rates in a thin-disk approximation, we hnd 
Cs oc ~ ,..-3/13^ depending on the optical depth of the disk and opacity sonrces (Paper 

I). This relatively flat temperatnre prohle results from the fact that for S oc r~^ (as in 
the onter parts) the potential is logarithmic and thus has a weak radial dependence. At 
a; 1, in contrast, this approximation may break down, since potential has stronger radial 
dependence. The isothermal approximation may not be justihed at the innermost region at 
later times (r -C Cgt). 



9 


A self-gravitating disk is locally stable, if 




CgK, 

ttGE 


> 


1 . 


(26) 


as long as the effects of viscosity and radial mass inflow are ignored (Toomre 1964).Here, k 
is epicyclic frequency and k, = 2^^‘^VL for H oc R~^. If we simply apply this criterion to the 
present model, we hnd Q ~ j at a; ^ 1 (Eq. 5 and 21), indicating that the disk 

is stable for HjR ^ 1/8. If HjR is small, gravitational instability will set out, making disk 
turbulent, thickening the disk (Paczyhski 1978). However, this is a very naive picture and 
a more sophisticated stability analysis, similar to Christodoulou et al. (1995a, 1995b) but 
including the effects of disk viscosity and radial gas inflow, is needed as future work. 


3.2. Formation of a Quasar Black Hole 


When Mj. exceeds a critical value at some radius, 

M„it(r) = (r/lO^-"cm)2/3M0, (28) 

the cloud will start to collapse due to a general relativistic instability (Shapiro & Teukolsky 
1983), resulting in the formation of a black hole. Equation (28) gives a critical mass (for a 
given radius) for spherical supermassive stars, while we are now concerned with evolution 
of a rotation-supported disk. Nevertheless, we employ the argument concerning spherical 
stars in order to see qualitative effects of general relativity, since a thin-disk approximation 
breaks down anyway near the center as mentioned above, and since a solid analysis of a 
collapsing self-gravitating disk based on the general relativistic formulation is not available 
at this moment. 

With this being kept in mind, we discuss a fate of a rotationally supported, viscous disk 
with a mass of ~ IO^Mq, a temperature of ~ 10"^K, and a size of several pc. In the present 
picture, such a relatively high disk temperature is preferable, since otherwise the disk will 
stay molecular rather than ionized. The accretion timescale is inversely proportional to the 
temperature, and hence it may exceed the age of the Universe for a molecular disk with 
a < 0.01 (e.g. Eq. 1 in Sasaki & Umemura 1996), unless alternative mechanisms, such as 
gravitational torque, remove the disk angular momentum. There are several possibilities 
to heat the disk. First, if the formation of primordial hydrogen molecules proceeds more 
slowly than the dynamical collapse, gas will not cool below ~ 10"^K. This may occur if 
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residual free electrons recombine quickly due to density enhancement, thereby suppressing 
the formation of a sufficient amount of H~ ions, which help to make hydrogen molecules 
(see Hutchins 1976, Palla, Salpeter, & Stabler 1983). Second, if the Universe was reionized 
through hrst-generation stars or objects, the disk will be effectively heated by strong UV 
background radiation (e.g. Sasaki & Umemura 1996). Finally, if star formation occurs 
within the disk itself, the disk material can be photoionized by stars. 

Figure 3 depicts the evolution of such a disk (by the solid lines) and the critical line 
for a gravitation instability (by the dotted line) in the (logr-logM,.) diagram. As time goes 
on, the disk becomes more and more condensed at the center, thereby increasing its mass 
within a hxed radius. The mass prohle is oc (Eq. 25) according to the self-similar 
solution, while the critical value gives Merit oc (Eq. 28). The solid line should cross the 
dotted line at 

Terit 0^ 10^^'® cm, Mr(rerit) OS lO'^'^ ^ M©. (29) 

We get a condensation of ~ lO^M© on a timescale of rsj 0.1 /-Q loV- 

The estimates above are optimistic, however, since it takes ro/cg lO^r to 
rQ/{a'cs) = 10®(Q;'/0.1)“^yr for accreting gas to reach the center, and thereby establishing 
a self-similar evolution of the disk. We thus safely conclude that within a timescale of 
~ 10^(Q;')“^yr a central region with a mass of lO^^^M© could become unstable, which may 
give rise to a proto-quasar black hole at high redshifts. Again, a general relativistic study 
of a collapsing rotating disk is necessary to conclude whether this scenario can work or not. 


We thank the referee for valuable comments and T. Hanawa and T. Tsuribe for 
useful conversation. This work is partial fulhllment of the Japan-US cooperative research 
program which is supported by Japan Science Promoting Foundation and National Science 
Foundation on US side. Analyses were, in part, made at Center for Computational Physics 
in University of Tsukuba, and Princeton University Observatory. Also, this work was 
supported in part by the Grants-in Aid of the Ministry of Education, Science, Sports, and 
Culture of Japan, 06233101, 08640329 (S.M.) and 06640346 (M.U.). 


A. Self-gravity under a thin-disk approximation 
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The most straightforward expression for the potentials under thin-disk approximation 


IS 


^(r) = 2G dd 

JO Jo 




T.{R)RdR 


(r^ + R^ — 2rR cos OyR ’ 


(Al) 


(e.g. Mestel 1963), where tq denotes the size of the disk and we ignored vertical mass 
distribution in the disk. After some algebra, we have 


— G{Ii + I2 + I3), (^ 2 ) 

where R, R, and R represent the Keplerian term, hnite contributions from the mass within 
R, and the mass beyond R, respectively, and are 


R = — [ 27iRJ]{R)dR, 
Jo 


R 


2n a 2 k 

k=l 


' {2k + l) 

Rk+2 Jq 


Jo 


4 = 27 r f; 021= (s(r) - 2fcr“-‘ R . 

with 

1 

a2k = - P2k{cos9)d9 = 

TT Jo 

{P 2 k is the Legendre function; see Eq. 24 of Mestel 1963). When S(r) 
particular, we hnd 

dT 27rGSoro 
dr r 


1 + ^ a2k 


fc=i 


2k 


(2fc)! 

{2>^k\y 


(A3) 

(A4) 


Soro/r, in 


(A5) 


We, hence, understand that if S(r) prohle is steeper than 1/r we may approximate 
gravitational attraction force to be —GiVA/r^ except near the outer edge. 
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Fig. 1.— Radial profiles of the self-similar variable, (t{x), as functions of a; = r/cst. The three 
solid lines represent the calculated values of a/a^ for a' = 10“^, 10“^, and 10“^, respectively, 
where Uo = a{x = 1). The transition radius at a; ~ a' separates the outer part, where 
a oc x~^, and the inner part, where a oc in each curve. Note that the dotted line 

corresponds to a/ao = x~^. 


Fig. 2.— Time evolution of a self-gravitating disk. /^From the top to the bottom, time 
development of M distribution, S prohle, and radial distributions of V^p (by the dashed line) 
and Vr (by the solid line). The units are tq = Ipc, Mq 2 Mgyr-i, Sq 30g cm 
Cs ~ 10km s“^, and to ~ lO^yr, respectively. Parameters are a' = 0.1 and ao = 50. The 
elapsed times are t/tg = 0.1 (indicated by i), 1.0, 10, and 10^ (indicated by f), respectively. 


Fig. 3. — Evolution of mass prohles of a self-gravitating disk with a total mass of 10®Mq, a 
temperature of lO^K, and a size of Ipc (by the solid lines). The attached numbers represent 
the elapsed times; t/to = 0.1,1.0, and 10. We assumed a' = 0.1 and ag = 50. Also displayed 
are the critical line for a general relativistic instability, rent (by the dotted line), and the 
Schwarzschild radius rg (by the short-dashed line), respectively. 



